Goal:

Think, explore, & write about what the co-evolutionary interaction between newts & snakes with different genetic architectures (GAs, combination of mutation rate & mutation effect size) can lead to. This markdown is investigation what is up with the different levels of correlation between rectangles and squares (in connection with GA1 tall). After fixing the row vs column error I looked at the correlation data and found that there was less correlation. So I decided to investigate why that might be and run a few more experiments. I am running an experiment to test how changing the square size might impact the calculations. I also plan on changing the interaction rate (but want to look at the math/ feasibility of it). This file contains results discussed in Tall_GA1!

Questions:

How does grid size impact the spatial correlation of newt and snake phenotypes?

Background

Experiment

I created a simulation study to observe the co-evolutionary outcome of the newt-snake interaction with different genetic architectures (GAs) in a spatial setting. I hypothesized that we would see an interaction (co-evolutionary arms race) between newt and snake phenotype under some GA combinations when newts and snakes were evolving over geographical space. Each GA is paired with another GA creating 16 combinations.

GA1 experiment values:

Landscape: A tall map!: 35*4 H, 35 W

I tested grid sizes for my tall landscape:

Each GA combination and trial has its own msprime simulation, but the msprime file is shared between the 5 section and 7 section slim runs.

The data

Preping the data by making singular files

## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## head: : No such file or directory
## All cor, lit, and grid files exist!
## This program will now end!

Loading in/ Combining the files to compare different grid sizes

generations

Correlation Histograms

In order to understand how spatial correlations where changing with time I took 5,000 generation time slices to look at all four trials correlation values. Each color is a different trial per GA combination. The histogram values are stacked. This section only looks at the 5 sections.

Plot 1

Plot 2

Plot 3

Plot 4

Plot 5

Plot 6

Plot 7

Plot 8

Plot 9

Plot 10

Phenotype Correlation across time

Next, we will examine three randomly chosen plots from both the 5 section and 7 section experiment. Time (in generations) in on the x-axis and both mean phenotype and phenotype spatial correlation in on the y-axis. Newt whole population mean phenotype is red, while snake mean phenotype is blue. The pink line is the phenotype spatial correlation.

Random 1 (5 sec)

Random 2 (5 sec)

Random 3 (5 sec)

Population Size Correlation across time

Random 1

Random 2

Random 3

Grid Figs

2

## [1] -0.6939078

3

## [1] -0.6396418

4

## [1] -0.5880404

5

## [1] -0.4995942

6

## [1] -0.4830822

7

## [1] -0.4421792

2 Pop Size

## [1] 0.703935

3 Pop Size

## [1] 0.8118469

4 Pop Size

## [1] 0.7261178

5 Pop Size

## [1] 0.5621264

6 Pop Size

## [1] 0.3736308

7 Pop Size

## [1] 0.3896028

## {-}